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I describe how the dynamics of galactic disks can be inferred by imaging and spectroscopy. Next I demonstrate that the 
decomposition of the rotation curves of spiral galaxies into the contributions by the various components of the galaxies 
is highly degenerate. Constraints on the decomposition can be found by considering implications for the dynamics of the 
galactic disks. An important diagnostic is the Toomre Q stability parameter which controls the stability of a galactic disk 
against local Jeans collapse. I also show how the density wave theory of galactic spiral arms can be employed to constrain 
the mass of a galactic disk. Applying both diagnostics to the example of NGC 2985 and discussing also the implied mass- 
to-light ratio I demonstrate that the inner parts of the galaxy, where the optical disk resides, are dominated by baryons. 
When I apply this method to the disks of low surface brightness galaxies, I find unexpectedly high mass-to light ratios. 
These could be explained by population synthesis models which assume a bottom heavy initial mass function similar to 
the recently proposed 'integrated galactic initial mass function' . 
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1 Introduction 

The dynamics of the disks of spiral galaxies is inferred from 
observations of the kinematics of the stars and the interstel- 
lar gas in the galactic disks. The most immediate way is to 
measure the mean rotation of the stars and the gas around 
the galactic center. Since the pioneering work of Rubin et 
al. (1978) rotation curves of galaxies have been observed by 
long slit spectroscopy in optical filter bands in great num- 
bers. The rotation curves are usually derived using emission 
lines which are emitted by HII regions in the spiral arms 
of the galaxies. In Fig. [T] I reproduce a part of a long slit 
spectrum of NGC 6070 close to the H/3 line (Fried & Fuchs, 
unpublished). The rotation curve - not yet corrected for in- 
clination - is clearly visible to the naked eye. Optical rota- 
tion curves reach radially outwards typically to three radial 
exponential scale lengths of the disks. 

Rotation curves observed by 21 cm emission from the 
interstellar HI, on the other hand, can usually be traced ra- 
dially much further outwards. In their pioneering work on 
NGC 3198 van Albada et al. (1985) could observe the HI ro- 
tation curve out to nearly 12 optical radial scale lengths and 
confirmed the existence of dark matter halos around spiral 
galaxies beyond any reasonable doubts. Rotation curves of 
both kinds are now available in large numbers and are com- 
piled in catalogues like that of Prugniel et al. (1998). 

In contrast to the mean velocities velocity dispersions 
are technically much more difficult to measure. The ob- 
served stellar absorption lines have to be deconvolved with 
the corresponding lines of template stars to measure line 
broadenings due to velocity dispersions of a few tens of 
km/s. Moreover spectra at at least two slit positions along 



Fig. 1 Long slit spectrum of NGC 6070 near the H/3 emis- 
sion line. The dispersion of the spectrograph is oriented hor- 
izontally, while the slit direction points vertically. 



the major and minor axes are required to disentangle the 
three components of the velocity dispersions in the radial, 
tangential and vertical directions, respectively. This disen- 
tanglement is greatly facilitated by the fact that the ratio 
of radial to tangential velocity dispersion can be related by 
epicyclic orbit theory (Lindblad 1959) to the shape of the 
rotation curves and need not to be measured directl}Q. Early 
results were presented in the work of Bottema (1993). How- 
ever, given the technical possibilities of that time his re- 
sults have too large errors to be used for the purposes de- 
scribed below. More reliable data have been obtained by 
Gerssen et al. (1997, 2000) for the two galaxies NGC 488 
and NGC 2985. Recently the disk mass project (Verheijen 
et al. 2004) has begun to measure the velocity dispersions 
of face-on galaxies with the accurate Pmass instrument. 
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2 Modelling rotation curves 

Rotation curves are modelled by constructing mass mod- 
els of the observed galaxies. These include usually a disk 
with an exponential density profile characterized by the cen- 
tral surface density and a radial scale length, a bulge model 
of the de Vaucouleurs or Sersic type characterized by the 
central density and the half-light-radius, and a model for 
the dark halo either with a NFW density distribution profile 
(Navarro etal. 1997), 
Po 




r(rc + r)^ 
or a quasi-isothermal sphere 

P(0 = o , o - 



(1) 



(2) 
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Fig. 2 Maximum disk decomposition of the rotation curve 
of NGC 3198. The contributions due to the stellar disk, the 
interstellar gas, and the dark halo are indicated. 
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As can be seen form eqns.[T]and[2lthe NFW law has a cen- 
tral density cusp, whereas the quasi-isothermal sphere has 
a central homogeneous core. Assuming essentially circular 
orbits of the stars and the interstellar gas around the galac- 
tic center - which is actually inadequate close to the center 
- the centripetal acceleration due to gravity is balanced by 
the centrifugal acceleration. For spherical systems like the 
bulge or the dark halo this simply means 

vl ^ GM(<t) ^ 

with G denoting the constant of gravitation, whereas for a 
flat exponential disk the circular velocity is given by Free- 
man's (1970) formula. The actual rotation curve of the galaxy| Fig. 3 Submaximal disk decomposition of the rotation 
is then given by the sum of three terms 




15 
R [kpc] 



(4) 



However, the decomposition of a rotation curve into its 
three components is notoriously degenerate. I demonstrate 
this again with the famous example of NGC 3198 which 
has a radial exponential scale length of rd = 2.6 kpc. In 
Fig. [lithe rotation curve is fitted by a maximum disk model, 
i.e. the disk contribution is chosen as high as allowed by the 
data. Fig. [3] shows a model with a submaximal disk, now 
of course with different halo parameters. Both fits cannot 
be distinguished by goodness-of-fit tests from each other, 
but further constraints on the mass of the disk are needed. I 
have suggested for some time (Fuchs 1999) that the implied 
internal dynamical state of the stellar disks according to the 
galaxy models might provide such constraints. 

3 Dynamical stability of galactic disks 

In a rotating self gravitating disk of stars there is an inter- 
play of gravity, pressure, and rotation. Patches of the disk 
become Jeans unstable on scales larger than the Jeans length 



curve of NGC 3198. 



angular momentum sets, on the other hand, an upper limit 
on the scales on which density perturbations can grow, 

(6) 



A < A. = 



where denotes the angular frequency of the patch under 
consideration. If the 'upper' limit is lower than the 'lower' 
limit, Ac < A J , growth of density perturbations is suppressed^ 
on all scales. More precisely this means that the Toomre 
(1964) stability parameter. 



Q = 



(7) 



3.36GS' 

must be larger than 1 to ensure that a thin galactic disk is 
stable against exponentially growing density perturbations. 
In eq. [7] denotes the epicyclic frequency of the stellar or- 
bits. 



V2^ 
r 



1 + — 



(8) 



A> Aj 



(5) 



where a denotes the velocity dispersion of the stars and S 
is the surface density of the disk. Detailed conservation of 



Vc dr ' 

and au the radial velocity dispersions of the stars. The coef- 
ficient 3.36 applies to a Schwarzschild velocity distribution 
of the stars; in the case of an isothermal gas it is replaced 
by TT. If Q < 1, this represents a serious threat to the fate 
of a galactic disk. In one of earliest numerical simulations 
of the evolution of a self-gravitating rotating disk Hockney 
& Hohl (1969) could already show that a thin Q < 1 disk 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2006) 



791 



stars gas stars gas 




— =10kpc 



Fig. 4 Snapshots of a numerical simulation of the evolu- 
tion of a combined disk of stars and gas. The disks are set 
up initially in radial force equilibrium and with Toomre pa- 
rameters Q* = 1.5 and = 0.5 — 0.7 respectively. Time 
is given in units of 10^ yrs. 

undergoes, even if it is rotating fast enough that the gravi- 
tational and centrifugal forces are in equilibrium, so violent 
instabilities that it eventually disintegrates. A Q = 1 disk, 
on the other hand, develops strong non-axisymmetric struc- 
tures, but stays essentially intact. More recently Fuchs & 
V. Linden (1998) have run numerical simulations of the evo- 
lution of a combined star and gas disk. Snapshots of one of 
these simulations are shown in Fig.lH Initially the combined 
disk was set up in radial force equilibrium and with Toomre 
parameters = 1.5 and Qg = 0.5 - 0.7. Such a disk has 
an effective Q parameter of slightly less than 1. The disk 
became immediately unstable as witnessed by the ring like 
structures especially in the gas disk. The rings fragment into 
clumps and shear then due to the differential rotation into 
spiral arms. These grow then rapidly to large amplitudes, 
both in the stellar and gas disks. The perturbations become 
so strong that they heat up the stellar disk dynamically so 
much that all spiral structures are eventually wiped out in it. 
A disk with Q < 1 is not an equilibrium model! 

The physical reason that galactic disks are prone to such 
fierce gravitational instabilities lies in their thin geometry. 
One can imagine a disk as an assembly of thin rings. The 
gravitational attraction does not vanish inside the ring like 
in a spherical shell, but diverges to infinite repulsion at the 
inner edge. At the outside edge of the ring the gravitational 
attraction becomes infinite regardless of the actual surface 
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Fig. 5 The radial gravitational force of a mass loaded thin 
ring calculated in the plane of the ring. The forces diverge 
at the edges of the ring. 




Fig. 6 Amplification factor of the swing-amplification 
mechanism for spiral density waves. The amplification fac- 
tor is shown as function of the circumferential wave number 
ky and the radial wave number kx, both measured in units 
of /ccrit, assuming a flat rotation curve, A = Q/2. Peak am- 
plification occurs Sitky = O.b /ccnt- 

density of the ring. In Fig. \5\ I illustrate an example. An 
assembly of such mass loaded rings must be obviously a 
highly unstable system. 

4 Density wave theory of spiral arms. 

The density wave theory predicts the preferential wave length] 
A of spiral density waves. I adopt here the concept of swing- 
amplification theory (Toomre 1981) which predicts 

A = X{A/n) Acrit (9) 
with the natural length unit 



The coefficient X in eq. [9] depends on the slope of the ro- 
tation curve measured by Oort's constant ^ = — ^ ^ ^ - In 
Fig.[6]l reproduce a figure of the operation characteristics of 
the swing-amplification mechanism taken from Fuchs (2001).| 
The amplification factor is shown as function of the cir- 
cumferential wave number ky = 2 tt/A and the radial wave 
number k^ assuming a flat rotation curve, A = Q/2. As can 
be seen amplification becomes maximal at k^ around 2 and 
3 /Ccrit( ^crit = 2 tt/ Acrit)- Peak amplification occurs in this 
example of a flat rotation curve at /Cy = 0.5 or A = 2 Acrit- A 
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S^crit 




Fig. 7 Sketch of how the preferred wave length - 2 A crit 
in this example - is related to the number of spiral arms. 

relation for the coefficient ?i{A/Q) for other slopes of the 
rotation curve is given in Fuchs (2001). 

Athanassoula et al. (1987) have pointed out that eq. [9] 
can be used to constrain the mass of the disk of a spiral 
galaxy. In Fig. [7] I illustrate with a sketch how the circum- 
ferential wave length is related to the number of spiral arms 
m. Obviously the number of arms is determined by how of- 
ten the preferred wave length fits onto the annulus, 



Thus one can infer from the number of spiral arms m, if the 
rotation curve is known, the critical wave length Acrit , and 
from thence the surface density H (cf. eq.[TO|). This concept 
has been tested successfully with numerical simulations of 
the spiral structure of galactic disks. Since the parameters 
of the disk which have been adopted for the simulation are 
known, one can read eq. [TT] from right to left and compare 
the predicted and actually found number of spiral arms. For 
instance, in the simulation by Fuchs & v. Linden (1998) 
the predicted number of spiral arms is m = 3 which can 
be clearly confirmed in the snapshot at t = 14.4-10^ yrs in 

Fig. a 

Often it is argued that spiral density waves do not grow 
spontaneously in the disks of galaxies, but are induced by 
tidal forces during encounters of two galaxies. Such events 
do occur, but are not long lived. Zang & Toomre (Toomre 
1981) have demonstrated with their work on the Mestel disk 
that such induced spiral density waves wrap up and die out 
on the timescale of one galactic revolution. Similarly Du- 
binski et al. (2008) have shown that in their numerical sim- 
ulations satellite galaxies orbiting around galaxies induce 
in the galactic disks of the parent galaxies transient, swing 
amplified density waves. 

5 Example: Diagnostics of NGC 2985 

NGC 2985 is a nearby (D = 1 8 Mpc) bright spiral galaxy. An 
image taken from DSS, which is reproduced in Fig.O shows 
that it has developed well defined bisymmetric spiral arms. 
Gerssen et al. (2000) have measured its rotation curve and 
the velocity dispersions of the stars in two slit positions so 
that the three dimensional velocity dispersion components 
are available. Moreover, they have imaged the galaxy and 
derived its surface density profile which is shown in Fig.O 



Fig. 8 Image of NGC 2985. (Source: DSS) 
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Fig. 9 Surface density profile of NGC 2985. The dashed 
lines are fits of the bulge and the exponential disk, respec- 
tively. The solid line is the superposition of both compo- 
nents. 

The authors have also provided a bulge-disk model for the 
photometry by fitting a bulge with an exponential density 
profile and an exponential disk with a radial scale length of 
Td = 2.6 kpc to their data. Arifyanto & Fuchs (2001) have 
changed this model slightly. They modelled the bulge by a 
softened power law, p ~ (r^ ^ ^2^-1.75 ^^^^ _ 275 
added a dark halo component and calculated rotation curves. 
In Fig. [Toll show first a maximum disk model in which the 
scale radii and the bulge-to-disk ratio as found in the pho- 
tometric model are adopted, but the mass-to-light ratio is 
chosen as high as allowed by the kinematical data. As can 
be seen in Fig. [Tol hardly any dark matter is required to fit 
the observed rotation curve. In the lower panel I show the 
predicted number of spiral arms and the Q parameter. Both 
diagnostics are consistent with the fact that the dynamical 
state of the disk of NGC 2985 allows the galaxy to develop a 
two-armed spiral. The implied mass-to-light ratio is M/ L b 
= 2.0 Mq/Lq. This is in good agreement with the predic- 
tion of population synthesis models. NGC 2985 has slB — R 
colour of B — R = lA mag. For such a colour the population 
synthesis models compiled by Bell & de Jong (2001) pre- 
dict a mass-to-light ratio of about M/Lb = 2.8 Mq/Lq. 

In Fig.[TT]l show a submaximal disk model for the rota- 
tion curve of NGC 2985. The fit to the rotation curve, which 
includes now a dark matter halo contributing as much as the 
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Fig. 10 Upper panel: Maximum disk model of the rota- 
tion curve of NGC2985. The contributions by the various 
components are indicated as dashed lines. Lower panel: Pre- 
dicted number of spiral arms m and the Toomre Q parame- 
ter. On the abscissa the scale is also indicated in arcminutes; 
1.5 arcmin corresponds to the optical radius of the disk. 
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Fig. 11 Same as Fig. [TOl but for a submaximal disk com- 
ponent. 



disk, is of the same quality as in the case of a maximum 
disk, but both the predicted number of spiral arms and the 
Q parameters as well as the implied mass-to-light ratio in- 
dicate that this not viable model for the disk of NGC 2985. 



6 Example: Diagnostics of the LSB galaxy 
F568-1 

In Fuchs (2002, 2003) I have applied similar diagnostics as 
in the previous section to a number of low surface brightness 
galaxies and describe here a typical result for the example of 
F568-1. The image of the galaxy, reproduced from de Blok 
et al. (1995) in Fig. [121 shows a regular two-armed spiral 
pattern. The rotation curve of F568-1 has been observed by 
de Blok et al. (2001). Unfortunately the velocity dispersions 




Fig. 12 Image of F568-1. 



of stars in LSB disks have not been measured up to now. De 
Blok et al. (2001) have provided - as for each galaxy in 
their sample - four types of models for the rotation curve 
of F568-1. The first model considers only a dark halo, the 
next includes the contribution by the interstellar gas, and 
the third model adds the contribution by a stellar disk with 
a, according to population synthesis, 'reasonable' mass-to- 
light ratio of M/Lr =1.4 Mq/Lq, which turns out to be 
a submaximal disk model. The fourth model is a maximum 
disk model. This suite of models is presented for the two 
cases that the dark halo is either a quasi-isothermal sphere 
or has a NFW density profile. Within each suite the reduced 
X^s are practically indistinguishable from each other. The 
surprising result is, when I constrain the decomposition of 
the rotation curve dynamically using the density wave argu- 
ment discussed in section (4), I find a mass-to-light ratio of 
M/L^=14! 



Lee et al.(2004) have attempted to modify the popula- 
tion synthesis model for LSB galaxies toaccommodate such 
high mass-to-light ratios. They treated the stellar popula- 
tion as an aging star burst with a bottom heavy IMF ^ ~ 
^-(i+a) (^^j^j^ ii^Q. notation that a = 1.8 corresponds to a 
Salpeter law). Moreover, they varied the metallicities of the 
stars in the range from [Fe/H] = -2.5 to 0.4. They find that 
they can indeed reproduce the mass-to-light ratios and the 
{B — R)o colours of about 0.8 mag, which are typical for 
the LSB galaxies in the sample of Fuchs (2003), by adopt- 
ing ages of the disks of 3 to 5 Gyrs, metallicities of [Fe/H] 
= -1.5 to -1, and an index of the IMF as steep as a = 2.35 
to 2.85, respectively. These parameters were chosen rather 
ad hoc to explain the results of Fuchs (2003). However, re- 
cently Kroupa and collaborators have developed the con- 
cept of an 'integrated galactic initial mass function' which 
predicts very steep IMF indices a. For instance, Weidner 
& Kroupa (2005) predict for galactic disks with masses of 
several 10^^ M©, which I find typically for the LSBs in my 
sample, in their 'maximal' model an a = 2.6 (in the nota- 
tion of Lee et al. 2004). The IGMF concept has not yet been 
adapted in detail to LSBs, but looks very promising and we 
hope to address this in the near future. 
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7 Conclusions 

Imaging and spectroscopy of galactic disks allows to infer 
their dynamics. However, the decomposition of the rotation 
curves of spiral galaxies is highly degenerate and needs fur- 
ther constraints. These can be found by considering the im- 
plications for the dynamics of the disks of the galaxies. A 
useful diagnostic of the dynamics of a galactic disk is the 
Toomre Q stability parameter. Another diagnostic is pro- 
vided by the density wave theory of galactic spiral arms 
which predicts the number of spiral arms. Both diagnostics 
can be used to constrain the masses of galactic disks. As an 
example of a high surface brightness spiral I have analyzed 
NGC 2985 along these lines and found that it seems to be 
dominated in the inner part, where the optical disk resides, 
by baryons. The resulting mass-to-light ratio is consistent 
with population synthesis models for such galaxies. The 
dynamical analysis of the disks of low surface brightness 
galaxies, on the other hand, implies unexpected high mass- 
to-light ratios. Such mass-to-light ratios can be explained by 
population synthesis models with a bottom heavy IMF. 
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